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ABSTRACT 

High-resolution N-body simulations of dark matter halos indicate that the Milky Way contains numerous 
subhalos. When a dark matter subhalo passes in front of a star, the light from that star will be deflected by 
gravitational lensing, leading to a small change in the star's apparent position. This astrometric microlensing 
signal depends on the inner density profile of the subhalo and can be greater than a few microarcseconds 
for an intermediate-mass subhalo (Mvh ^ 10'* M©) passing within arcseconds of a star. Current and near- 
future instruments could detect this signal, and we evaluate the Space Interferometry Mission's (SIM's), Gaia's, 
and ground-based telescopes' potential as subhalo detectors. We develop a general formalism to calculate a 
subhalo's astrometric lensing cross section over a wide range of masses and density profiles, and we calculate 
the lensing event rate by extrapolating the subhalo mass function predicted by simulations down to the subhalo 
masses potentially detectable with this technique. We find that, although the detectable event rates are predicted 
to be low on the basis of current simulations, lensing events may be observed if the central regions of dark matter 
subhalos are more dense than current models predict (> 1 Mq within 0. 1 pc of the subhalo center). Furthermore, 
targeted astrometric observations can be used to confirm the presence of a nearby subhalo detected by gamma- 
ray emission. We show that, for sufficiently steep density profiles, ground-based adaptive optics astrometric 
techniques could be capable of detecting intermediate-mass subhalos at distances of hundreds of parsecs, while 
SIM could detect smaller and more distant subhalos. 
Subject headings: astrometry - dark matter - Galaxy:halo - gravitational lensing:micro 



1. INTRODUCTION 

Numerical simulations of dark matter halos have re- 
vealed the presence of numerous subhalos over a wide 
range of masses ex tending down to the simulations' res- 
olution limits ( e.g., iGhigna et alj [19981 lIOOOl: iKlvpin et al 



1999a; Moore etail 119991: iKravtsov et all I2004t iGao et al 



2004; Reed et al."2 005l; iDiemand et alj|2007bl; iSpringel et al 



2008; Diemand et all l2008h . This substructure is the rem 
nant of hierarchical structure formation; as halos merge to 
form larger structures, the inner portions of the ancestor halos 
become subhalos. If all halos leave subhalo remnants, then 
the subhalo mass function may extend to ma sses far smaller 
than the resoluti on limit of any simu lation (iHu et al J 120001 : 
IChen et al.ll200HlProfumo et al.ll2006l:lbiemand et alJl2005h . 
Although subhalos may be destroyed by gravitational in- 
teractions with the host halo, other subhalos, and stars, 
there are indicat ions that their dens e inner regions sur vive to 
the present dav (jHayashi et alJl2003 : Kazantzidis et al Jl2004 



tne present aav (|Jriayasm et aiJlzuuj ; ivazantziais et ai jizuu4t 
van den Bosch et al."2005';'Read et al. 2006; Berezinskv et alJ 
2006; Zhao et al. 2007; Green & Goodwin 2007; Goerdt et al] 
2007t [Schneider et alj 1201 ft llshivama et alj 12010). More- 



over, high-resolution simulations of halos similar to the Milky 
Way's host suggest that subhalos with masses greater than 
10^ Mp) ai'e pre sent at the Solar radius (ISpringel et al J 120081: 
IDiemand et al.l |2008). Although simulations indicate that 
these large subhalos are significantly disrupted by the Galac- 
tic disk (D'Onghia et al. 2010; Romano-Diaz et al. 2010), the 
resolution is not sufficient to determine the fate of the sub- 

' Canadian Institute for Theoretical Astrophysics, University of 
Toronto, 60 St. George Street, Toronto, Ontario M5S 3H8, Canada; 
erickcek@cita.utoronto.ca 

^ Perimeter Institute for Theoretical Physics, 31 Caroline St. N, Water- 
loo, Ontario N2L2Y5, Canada 

^ Dunlap Fellow, Dunlap Institute for Astronomy and Astrophysics, 
University of Toronto, 50 St. George Street, Toronto ON M5S 3H4, 
Canada 



halos' innermost regions. It is therefore possible that numer- 
ous subhalos are located within a few kiloparsecs of the Solar 
System, with profound implic ations for both direct detection 
of th e dark matter particle (jKamionkowski & KoushiappasI 
120081) an d indirect detectio n through its annihilation signa- 
ture (e.g., Bergstrom et al.''T999; Calcaneo-Roldan & Moor3 
pOO; Stoehretal. 2003; Diemand et al. 2007a; Ando 2009F 
iKamionkowski et al. 20 IC^. 

Unfortunately, both indirect and direct detection of dark 
matter continues to be elusive, and gravity still provides 
the only uncontested evidence for dark matter Gravita- 
tional lensing is an especially powerful tool in the study of 
dark matter substructure; subhalos within our galaxy could 
be detected t hrough their effec ts on signals from millisec- 
ond pulsars (ISiegel et alJ 120071) . and substructure in lens- 
ing galaxies has several observational signatures. Subha- 
los have been proposed as the origin of observed flux-ratio 
anomalie s between th e multiple images of strongly lensed 
quasars dMao & Sch neider 1998; Metcalf & Madau 200^ 
Chiba 2002; Dalai & Kochanek 2002; Kochanek & Dal3 
2004)' they can a lso alter the time delays between these im- 
ages (IKeeton & M oustakas 2009; Congdonetal. 201( j) and 
their separations (Koopmans et al. 2002;; iChen et al.l 120071 : 
IWilliams et al. 2008; More et al. 2009). If flie source is 
extended, then subhalos ca n distort the image's shape an d 
surface brightnes s (Metcajg 120021; llnoue & Chib^ l2005allbl : 
iKoopmansI l2005t l^getti & Koopmans 2009a b). Intrigu- 
ingly, the observed lensing anomalies can only be explained 
if the central regions of the lensing galaxies contain signifi- 
cantly more substruct ure than is predicted by n-body simula- 
tions LM ao et al. 2004lAmara et al. 2006: Maccio e t al.ll2006l: 
IXu et al . 2009, 2010), although it has be en sugg ested that 
these studies use aty pical lensing galaxies (iBryan et al.ii2008l : 
iJackson et alllMTol) . 

Individual subhalos within a lensing galaxy could be de- 
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tected if they strongly lens one quas ar image, splitting it 
into a closely spaced pair of images jYonehara et alj 120031 : 
llnoue & Chiball2005al) . Unfortunately, the diffuse nature of 
dark matter subhalos implies that they have small Einstein 
radii (Zackrisson et al. 2008). Consequently, the split images 
are resolvable only for the largest subhalos, making it unlikely 
to find a large enough subhalo with a sufficientl y small im- 
pact p arameter to detectably split a quasar image tRiehm et alj 



pact p a 
I2OOI . 



In this paper, we consider a different way to find subhalos 
through gravitational lensing: instead of looking for split im- 
ages, we investigate how the astrometric deflection of an im- 
age changes as a subhalo moves. This astrometric microlens- 
ing approach has two advantages over strong lensing. First, 
it is much easier to measure a change in the position of the 
centroid of an image than it is to resolve an image pair into 
distinct sources. The minimal separation required to resolve 
a pair of images is limited to approximately the resolution of 
the telescope (> 25 mas), whereas at high signal-to-noise ratio 
(S/N), the position of the centroid of the image can be mea- 
sured with hundreds of times higher precision. Second, strong 
lensing only occurs when the angular separation between the 
source and the lens is smaller than the lens's Einstein angle, 
while astrometric deflection is detectable for far larger impact 
parameters. The disadvantage associated with astrometric mi- 
crolensing is that it must be a dynamical event because we do 
not know the true position of the source. As the lens moves 
relative to a background star, the star's image will move as 
well, and that is the detectable signature of astrometric mi- 
crolensing. We are therefore constrained to local subhalos 
with significant proper motions. 

Astrometric lensing signals from dark matter subhalos are 
necessarily small because subhalos are diffuse, and astrome- 
try has only very recently progressed to the point that an astro- 
metric dark matter search is feasible. The rapid development 
of 1-100 |/as astrometric precision has been driven by a wide 
variety of fie lds - for examp le, following orbits in the galactic 
center (e.g., iLu et alJlIOO^ [Gillessen et al. 2009); astromet- 
ric detecti on of pla nets (e.g., Unwin et al. 2008; iLaw et akl 
12009; MalbetetaL J2OIO); accurate parallax determination 
(e.g., Hen ry et alJ l2009t iSubasava ge et al. 200 9b; and flie 
determination of st el lar orbits ("e.g.. . Helminiak etaP 120091; 
Pravdoetal. 2005: iKonopackv et alj 120101; [Ireland et all 
2008uiDupuvetalJl2009h . 



New space-based a strometric missions such as Gaia 
(iLindegren et alJ l2008h and SIM (lUnwin et al.1 llOOSh are 
opening the possibility of ultra-high-precision all-sky and 
targeted searches. Targeted ground-based astrometry is al- 
ready capable of lOO/^as precision in arcminute-sized fields, 
and ne w larger telescopes wi ll significantly improve that pre- 
cision (Cameron et al.l l2009h . while even higher precisions 
are possible on bright stars (e.g., Muterspaugh et al. 2006; 
Ivan Belle et alJ l2008). 

These technical advances have inspired considerable inter- 
est in astrometric m i crolensing by stars and dark compact ob- 
iects (Walkei^' 19951; 'Paczvnski"1995', '1998'; iMiralda-Escud^ 
[1996; Gould 2000; Gaudi & Bloom 2005), and by baryonic 
clouds (Takaha shil 120031: ILee et al.ll2010l) . In this paper, we 
investigate if these instruments can also be used to search for 
dark matter subhalos. We find that, for standard dark matter 
models, observing a subhalo lensing event rate during a blind 
astrometric survey is highly unlikely. If the central regions 
of dark matter subhalos are denser than expected, however, 
the lensing event rate can be significantly enhanced. We also 



explore the possibility of using astrometric lensing to con- 
firm the presence of a subhalo detected through its gamma-ray 
emission. We find that ground-based telescopes could detect 
lensing by a nearby (~50 pc) subhalo with a post-stripping 
mass greater than 1000 M©, while SIM could probe these sub- 
halos at greater distances (~100 pc) and detect nearby subha- 
los with one-tenth of this mass. 

This article is organized as follows. In Section |2] we de- 
scribe the image motion induced by subhalo lenses with a va- 
riety of different density profiles, namely, a singular isother- 
mal profile, the NEW profile, and a generalized power-law 
density profile. In Section [3] we develop an astrometric ob- 
serving strategy designed to reliably detect subhalo lensing 
while rejecting possible false-positive detections. In Section 
m we calculate the areas of sky over which particular sub- 
halo lenses are detectable and determine the all-sky subhalo 
lensing event rates for several models of dark matter substruc- 
ture. In Section |5] we evaluate current and planned astromet- 
ric survey capabilities in the context of detecting lensing from 
a subhalo in both all-sky and targeted searches. We summa- 
rize our findings and conclude in Section |6l To evaluate the 
lensing signatures and event rates, we developed models of 
the concentration-mass relation and the mass function for lo- 
cal su bhalos based on the findings of the Aquarius simula- 
tions ( Springel et alj|2008l) : these models are presented in the 
appendix. 

2. ASTROMETRIC SIGNATURES OF LENSING BY SUBHALOS 

The shapes of the intermediate-mass subhalos that we will 
be considering cannot be probed by current numerical simu- 
lations; to simplify our analysis we will assume that the sub- 
halos are spherically symmetric. Numerical simulations can 
only probe the shap es of the larg est subhalos with masses 
greater than 10^ Mg (iKuhlen et alJl2007l: iKnebe et al.l l2010^. 
While these subhalos are triaxial, they tend to be more spher- 
ical than thei r host halos. In a simulation that includes bary- 
onic phvsics. TKnebe et al] (12010.) found that subhalos located 
within half of the virial radius of the host halo have nearly 
spherical mass distributions in their innermost resolved re- 
gions, with a median minor-to-major axis ratio of 0.83. We 
will see that astrometric microlensing is only sensitive to the 
density profile near the center of the subhalo (r < 0.1 pc); 
as long as the inner region of the subhalo is nearly spherically 
symmetric, the subhalo's triaxiality is irrelevant for our analy- 
sis, and we do not expect deviations from spherical symmetry 
to significantly affect our results. 

We will also assume that the dark matter subhalo is trans- 
parent and that it contains no stars. The presence of a star 
in the center of the subhalo would change its lensing signa- 
ture if the star's mass is comparable to the dark matter mass 
in the central 0.1 pc of the subhalo, but stars further from 
the center would not have an effect. Furthermore, the dis- 
crepancy between the number of intermediate-mass subha- 
los seen in simulations and the number of dwarf spheroidals 
observed i n the Milky Way, known as the " missing satellite 
problem" (jKlvpinet al.lll999bl: iMoore et al][l999,) . indicates 
that stars are rare in subhalos with masses less than 10^ M© 
(iMadau et alJ 1200 8). Finally, we will assume that the sub- 
halo's diameter is small compared to both the distance d]^ be- 
tween the lens and the observer and the distance c/ls between 
the lens and the star. 

Figure[T]shows a schematic view of lensing by an extended 
transparent object like a dark matter subhalo. When a light ray 
passes through a spherically symmetric thin lens, the image of 
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Fig. 1. — a) Diagram showing tiie position of tlie source star (in black), its image (in gray), and the lens (gray ellipse). We will generally assume that a <*: jS so 
that the ray's impact parameter in the lens plane (f = d^d) may be approximated as f ^ d^fi. b) The same lensing system viewed as projected on the sky. The 
center of the lens is moving with velocity vx along the x-axis. 



the star is shifted from its true position by an angle 



ds 



4GM2d(0 



^, 



(1) 



where ds is the distance between the observer and the star, ^ 
is the distance between the center of the lens and the star's 
image in the lens plane (^ = dj^O), and ^ = ^/^ points from 
the lens to the star Throughout this work, we set the speed of 
light c - I. The mass M2d(^) in Eq. ([TJ is the mass enclosed 
by the cylinder interior to ^ and is obtained by integrating the 
projected surface mass density S over the area of the circle 
with radius ^. 
As the subhalo moves relative to the background star, the 

angle /3 that extends from the lens to the star will change, 
and the position of the star's image will change accordingly. 
We take the star to be fixed at the origin of an xy coordi- 
nate system on the sky, and we define the x-axis to be parallel 
to the subhalo's transverse velocity vt, as shown in the right 

panel of Fig.[T| The vertical component of yS is therefore fixed 

l^_,,(f)=y6,,o],and 



y8.«=/3.,o-4.2" 



Vt 



200km/s 



50 pc 



t 
5^- 



(2) 



where /3^fl is the value of /3x at f = 0. We see that a nearby 
subhalo will move several arcseconds during a multi-year ob- 
servational period. 

For the subhalos we consider, the deflection angle a will 
be on the order of microarcseconds. Since /? changes by sev- 
eral arcseconds over the course of the observation, /3 '» a for 
most of the observational period. We will further assume that 
y6y,() » a so that we are always considering the weak-lensing 
regime with/? » a. We verify in Appendix lAl that this con- 
dition is satisfied for all subhalo lensing scenarios, provided 
that d^ <^ 1000 kpc. This confirms that we are firmly in the 
weak lensing regime as long as we only consider subhalos in 
our local group. In this case, there is only one image of the 
star, and it is always located on the line connecting the lens 
position to the star's position, with the star between the image 
and the lens. We will use the /3 » a assumption to simplify 
the lens equation by approximating p ^ 6. In this case, ^ may 
be approximated as ^ ^ diJ3, and Eq. ^ becomes a simple 
equation for the deflection angle a in terms of the impact pa- 
rameter y8. In the following subsections, we will use Eq. ([TJ to 
show how the path taken by the star's image during a subhalo 
transit depends on the subhalo's density profile. 



2.1. Singular isothermal sphere 
The density profile for a singular isothermal sphere (SIS) is 



pW = 



InGr^ 



(3) 



where a-y is the velocity dispersion of the halo. Although 
numerical simulations indicate that large dark matter subha- 
los without baryons do n ot have this steep a density profile 
dSpringel et al. 2008; Die mand et al J 120081) . we consider the 
SIS case in detail because it simply illustrates key features 
that are shared by the astrometric lensing signatures from dark 
matter halos with shallower profiles. 

The two-dimensional enclosed mass for an infinite SIS is 
^2d(^) = ^crl^/G. Since M2D depends linearly on ^, Eq. ([TJ 
reveals that a is independent of the separation between the 
lens and the star The deflection angle is always the Einstein 
angle of the SIS: 



r^ 



d] 
ds 



^) ^""'' 



lOyuas 



cr,, 



0.6km/s 



1 _ ^ 
ds 



(4) 
(5) 



There are two images, with d - +6^^0, only if jS < 6^^. 
We will only consider larger impact parameters, in which case 
there is only one image, with a - 6^^^. As the SIS moves rel- 
ative to the star, the direction of the deflection angle changes. 
For an infinite SIS moving from the distant left to the distant 
right, the image starts 6^^ to the right of the star's true posi- 
tion and then traces a half-circle with radius G^^ until it ends 
6^^ to the left of the star's true position. 

For an SIS, the mass enclosed in a sphere of radius R is 
proportional to R; if the SIS has infinite extent, then its mass 
is infinite. It is customary to characterize an SIS by its virial 
mass: the mass contained in a sphere with mean density equal 
to the virial density pvir- iBryan & Norman (1998) provide a 
fitting formula for the virial density in a flat ACDM universe, 

Pvir = {l8n^ + 82[Qm(z) - 1] - 39[Qm(z) - l]')pcrit, (6) 

f^M0(l +zf 



^m(z) 



(7) 



f^Mo(l +zP-h1-Qmo' 
Pcntiz) = (o.0924 ^ j h^ [Qmo( 1 +z)' + l- Omo] , (8) 

where Hq = lOO/ikms 'Mpc ' and Qmo is the present-day 
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matter density divided by the critical density. We will use 
standard cosmological parameters: h - 0.7 and Quo = 0.3. 
The virial density is 4.6 Mg pc"^ at redshift zero and it in- 
creases monotonically with redshift. The velocity dispersion 
cr^, in terms of the virial mass is 



oi.^G 



7rpvir(z,)M; 



2 \ 



1/3 



(9) 



where z,, is the redshift at which the halo's virial mass is eval- 
uated. To facilitate comparisons with N-body simulations, we 
take Zv = in our calculations, but we note that increasing 
Zv would make the subhalos denser and would enhance their 
lensing signals. Inserting this expression into Eq. (|4]i gives 



9f= (l.l;aas)(l-^ 



Mv 



2/3 



Pvir(Zv) 



1/3 



104 Mq/ \4.6Mepc- 

(10) 
These properties describe an SIS in isolation. Once a 
subhalo is accreted by a larger halo, the outer tails of its 
density profile are stripped of mass. Numerical simula- 
tions indicate that a subhalo in our Galaxy may lose be- 
tween 99% and 99.9% of its initial mass due t o tidal strip- 
ping fr om the smooth compone n t of the halo dHavashi et alj 
2003t iKazantzidis et alJ l2O04t Ivan den Bosch et al.l 120051 
Readetal.ll2006^ . and stars will further s t rip the outer por- 
tions of subhalos ('Berezinskv et al."2006'; 'Zhao et al."2007^. 
Green & Goodwin 2007; Goerdt et al. 2007; Schneider et al. 
2010t llshivama et al. 2010). We will deal with this trunca- 
tion by defining a truncation radius /?, and setting p = for 
R > R,. The mass contained within /?, is the mass of the sur- 
viving subhalo Mf. We will describe the tidal stripping with 
the parameter m\,i = M,/Mvii-, where Myir is the original virial 
mass of the subhalo, evaluated at redshift z,,. The truncation 
radius is then given by 



R, = (0.56 



^ ^\0.001/\ 



Mvi 



104 Mf 



1/3 



4.6 Me pc 



-3\l/3 



Pvir(Zi.) 



(11) 



The angular size of the truncated halo is 



0, = 0.64' 



/ 'Wbd \ 

1 0.001/ 



50 pc 



Mvi 



IO^Mf 



"''4.6M^pc-^^''' 



Pvir(Zi) 



Thus we see that 6^^ <& 6, for all subhalos of interest. 



(12) 



For a truncated SIS, the two-dimensional enclosed mass is 



M2d(^ < R,) 



2(tI 
G 



^tan 




^ 



■e 



(13) 

and M2D = M, for f > iJ,. If ^ > Rt, then the deflection angle 
is the same as for a point mass with Einstein angle 



^.^^ 



4GM,^LS 
dhds 



9f f 



(14) 



Thus we see that (^^ <K 6, implies that 0™ <K 0,. Therefore 
we may approximate the position of the brightest image as 
a = (0l^Y//3 for/? > 0,. When we insert Mid from Eq. O 
into Eq. dU for a and assume that ff <s; /?, we find that the 
deflection angle for lensing by a truncated SIS is 



,sis 



cJ^-e^yS 



f^(f)f- 



qSIS 



<j6< 



for yS > 



(15) 
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Fig. 2. — The deflection angle for a star at 5 kpc as an SIS subhalo 50 pc 
away moves at a velocity of 200 km/s from the far left to the far right. This 
subhalo's velocity dispersion is 0.72 km/s, which corresponds to a virial mass 
of 5x10 Mq, but it has been stripped to a radius of 0.02 pc and contains 
only 5 /Wq. In the top plot, the subhalo center passes 1 arcsecond below the 
star, and the circles show the image's position every year for the 100 years 
surrounding the time of closest approach. In the bottom plot, the subhalo 
center passes 50 arcseconds below the star, and the circles show the image's 
position every 10 years for the 300 years surrounding the time of closest 
approach. 
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Figure |2] shows how the image of a fixed star moves as a 
truncated SIS subhalo passes below the star's true position on 
the sky. The lens's Einstein angle is 0|'^ - 15/ias, which cor- 
responds to cTy - 0.72 km/s and Myir = 5x10^ M©. The image 
motion is highly sensitive to the ratio vx/iiL because this ratio 
determines how/3 changes with time [see Eq. (|2]l]. If vt/i^l 
is decreased, then the change in /? during a set time interval 
is decreased, and the image motion slows down. Through- 
out this work, we adopt vx = 200 km/s because that is the 
typical velocity of a dark matter p article in the halo (e.g., 
iDrukier et alJll986l; IXue et alj|2008l) . Wifli fliis velocity, we 
will see that a detectable subhalo must have d\^< 100 pc, and 
we adopt c/l = 50 pc as our fiducial lens distance. Provided 
that c/l *^ i^s, the distance to the source has a minimal im- 
pact on the image motion because lis only enters through the 
factor (1 - d]^lds) in 0^^ [see Eq. dTOlll. We use d^ - 5 kpc 
as our fiducial value because this is a reasonable distance to a 
target star; changing d^ to any value greater than 1 kpc would 
have no noticeable effect on the image motion. To illustrate 
the effects of subhalo truncation in Fig.|2l we assume that the 
subhalo has been extremely stripped by close encounters with 
stars so that its radius is 0.02 pc (0, - 85"), which implies 
that M, = 5 Mq. We see that the image motion consists of 
three stages: as the edge of the subhalo approaches from the 
far left, the image very slowly moves rightward until the star 
is behind the subhalo, then the image rapidly traces an arc 
as the subhalo center passes by the star, and finally the star 
slowly returns to its true position as the subhalo moves off to 
the right. 

The impact parameter j3y determines how quickly the im- 
age moves as the subhalo passes by the star. In the top half of 
Fig. |2]the lens impact parameter is 1", and the image rapidly 
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Fig. 3. — The deflection induced over 4 years by a moving SIS lens with 
the same properties as in Fig.|2] The path of the lens is depicted by a dotted 
arrow. To show the image trajectories, the image motion is exaggerated a 
factor of 10* relative to the star's positions; a scale bar corresponding to 20/ias 
motion is shown. Twenty equally spaced measurement points over the 4-yeai' 
observational period ai'e shown for each curve. Note that the stars closest to 
the lens position undergo much more rapid position changes. 

traces out a semi-circle of radius 0|'^ during the few years 
surrounding the time of closest approach, just as if the lens 
had infinite extent. The effect of the SIS's truncation is more 
apparent in the bottom half of Fig. |2l where yS, = 50"; the im- 
age's trajectory is closer to a circle, and it will become more 
and more circular as /?,, increases. The image transverses this 
circle very slowly, taking 10 years to move 2 //as, in contrast 
to the image in the top panel, which moves nearly 30 yuas in 
only 5 years. Thus we see that the only detectable portion of 
the image's path in the sky is the period surrounding the mo- 
ment of closest approach between the star and the lens, and a 
small impact parameter is required to make the image move 
significantly over the course of a few years. Figure [3] further 
illustrates the necessity of a small impact parameter by show- 
ing how the images of stars at different positions relative to 
the lens move over the course of 4 years; only the stars along 
the lens's path with/?,. < 2" have images that are significantly 
moved during the observational period. For stars that are this 
close to the center of the subhalo, with/? <s Q,, the truncation 
of the density profile does not affect the image trajectories, 
as seen in the top panel of Fig. |2] We will therefore assume 
thatjS <K Q, for all interesting lensing scenarios and ignore the 
subhalo's truncation when considering other density profiles. 

2.2. NFW density profile 
The NFW profile. 



P(r) 



m-tr 



(16) 



was found to be a universally good fit to the density profiles 
of galaxy and cluster halo s in early numerical simulations 
dNavarro et al.lll996l I1997D . The two-dimensional enclosed 
mass for a subhalo with virial mass Mvir and an NFW density 
profile with concentration c = Ryu/rs is 



M2D© = 



Mv, 



ln(l+c)-^ 
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(lBartelmann "T996; Golse & Kneib 2002). 

From Eq. ([T]), we see that the magnitude of the deflection 
angle a is proportional to M2d/^- For the NFW profile, a oc ^ 
if r <s: rs, and a oc ^ ' if r » r,. Therefore, as an NFW sub- 
halo approaches a star, the deflection angle will increase until 
the star crosses the scale radius (^ ^ r,), and then it will de- 
crease until the star crosses the subhalo center, at which point 
it will begin to increase again until the star crosses r, on the 
other side of the subhalo. In this sense, the scale radius of an 
NFW profile plays the same role as the truncation radius for 
a truncated SIS. If the impact parameter is close to the scale 
radius (J3y - rjdi), then the image trajectory is roughly circu- 
lar, and it resembles the bottom half of Fig. |2] Unfortunately, 
the subhalos that are massive enough to deflect the star's im- 
age by several microarcseconds (Mvh- <: 10'^ Mq) have large 
scale radii (r, > 2 pc for c < 100); /?,. - '".v/'^l is a large im- 
pact parameter, and the image position changes very slowly 
as the subhalo moves. Moreover, just as with a truncated SIS 
lens, the reversal in the image's motion as the star crosses the 
scale radius (J3 ^ r^jdi) is very slow, regardless of the impact 
parameter ySy 

As in the SIS case, the most promising lensing scenario oc- 
curs when the center of the NFW subhalo passes very close 
to the source. The key difference is that a is nearly constant 
for ^ <K 7?, if the lens is an SIS, which leads to the semi-circle 
image trajectory displayed in the top portion of Fig. |2l For 
an NFW lens with f <s rj, the deflection angle is very small, 
as shown in the bottom panel of Fig. |4] The NFW density 
profile leads to a no-win situation: if you decrease the impact 
parameter jSv in order to enhance the change in the image's po- 
sition over a set time period, the magnitude of the deflection 
decreases. We are forced to conclude that astrometric lensing 
by subhalos is only detectable if the inner density profile is 
steeper than p oc r"' . 

2.3. Generalized density profile 

We have seen that astrometric gravitational lensing by sub- 
halos is only detectable if the center of the subhalo passes 
close to the star's position during the observational period, 
during which the subhalo moves about 0.001 pc (for a 5-year 
observational period). Therefore, only the innermost por- 
tion of the subhalo is responsible for the astrometric lens- 
ing signature. Unfortunately, very little is known about 
the intermediate-mass (10 Mq < M, < 10'' Mq) subhalos 
that are capable of producing detectable astrometric lensing 
events. High-resolution N-body simulations can probe the 
density profiles of only the largest {M, > 10** Mq) subha- 
los, and even these profiles are unresolved at r < 350 pc 
(lSprinp:eleta l.ll2008; Dieman det al...2008.) . For these large 
subhalos, JDiemand et al. (200^ find that p oc r ~^-^ in the 
innermost resolved regions, while ISpringel et alj (12008!) see 
p oc ^-(i-2-i-7) 2ii their resolution limit for nine large subhalos, 
with no indication that the slope had reached a fixed central 
value. Meanwhile, at the opposite end of the mass spectrum, 
iDieman d et al. (2005) find that the first Earth-mass dark mat- 
ter microhalos have steeper density profiles with p oc ^-(15-2.0) 
at redshift z - 26, and higher-resolution simulations indicate 
that this steep profile extends to within 20 AU of the micro- 
halo center (Jshiyama et al.1,2 010,) . 
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Fig. 4. — The deflection angle generated by a moving lens with rfs = 5 kpc, 
^L = 50 pc and vx = 200 km/s. The virial mass of the lens is 5 X lO' M©, 
and its concentration is Svi,/r_2 = 99. The inner density profile of the lens 
is given byp oc r^~f ^ and the dift'erent panels correspond to dift'erent values of 
y. The impact parameter is 1 arcsecond, and only the portion of the image 
path con'esponding to the time surrounding the moment of closest approach 
between the image and lens is shown. Note that the image path becomes more 
linear and the image motion slows down considerably as y is decreased. 

In light of this uncertainty, we consider a generic density 
profile 

/ ,-\-y 

(18) 



P^X) = PC — 
\'-o 



with 1 < y < 2. We assume that a constant-density core, 
if present, is significantly smaller than our typical impact pa- 
rameters of 0.001 pc, and we assume that the subhalo does not 
contain a black hole. Larger cores would decrease the lens- 
ing signal while the presence of a black hole would enhance 
it by adding a point mass and steepening the density profile 
dBertone et'al. 2005; Ricotti & Gould 2009). If we take this 
density profile as infinite when calculating the projected sur- 
face density E, we find that 



2(^) = V7rpo'"o 



r[0.5(y-l)] /^ 
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where Y\x\ is the Euler gamma function. 

Of course, this density profile does not extend to infinity; 
the subhalo's density profile will be truncated by tidal strip- 
ping, and it may also transition to a steeper power law, as in 
the case of an NFW profile. If the density profile is truncated 
at r = R,, then the surface density diverges from Eq. ( fT9] l as £, 
approaches R,, but for^ <k /?,, Eqs. ( fT9] l and (|20] | are still good 
approximations. For instance, if y = L5(1.2), M2d(^) for a 
subhalo truncated at R, is greater than 80% (50%) the value 
given by Eq. ( |20] | if ^ < 0.1/?,. We will show in Appendix 
|A]that detectible astrometric signatures are only produced if 
£, < 0.03 pc, and Eq. ( |20| ) is accurate to within 20% for sub- 
halos with y > 1.5, Mvk < 10'^ Mq, and Ri > 0.1 pc. Fur- 
thermore, the lower bound on R, is significantly smaller for 
subhalos with Myir <s 10^ Mq. We will therefore use Eq. ( l20l i 
and take /?, > 0.1 pc as a conservative lower bound, although 
we note that the resulting deflections may be slightly overes- 



timated, especially if y < 1.2. As shown in Fig.H] however, 
detecting a subhalo with y < 1 .2 is challenging, and we con- 
clude that Eq. ( |20] | is accurate to within ~ 20% for subhalos 
of interest. 

If a dark matter subhalo with a density profile given by 
Eq. (fTsT l passes in front of a star, Eq. ([T]i tells us that 
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where we have defined 
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(22) 

Like 6^^, 6„ depends on the distances to the lens and the 
source only through the factor (1 - dL/ds). We also note that 
6a is related to the Einstein angle 9e'. 
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(23) 



We will continue to assume that a <^ /3 so that ^ (see Fig.[TJ 
is approximately equal to diJ3. 

Equation ( l22b gives the magnitude of the deflection angle 
in terms of the parameters of the density profile ro and po, but 
this is not the most useful description of the subhalo. Instead 
we characterize the subhalo by either its mass after tidal strip- 
ping (M, = njbdA^vir) or the mass contained within a radius 
of 0.1 pc from the subhalo center (Mo.ipc). Although M, is a 
more standard and intuitive description of the subhalo mass, 
using Mo.ipc offers two advantages. First, Mo.ipc completely 
determines the deflection angle; without loss of generality, we 
can set ro = 0.1 pc, in which case 
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(24) 

Second, Mo.ipc is the portion of the subhalo's mass that is ac- 
tually probed by astrometric microlensing because truncating 
the subhalo's density profile at R, - 0.1 pc does not affect 
its astrometric lensing signature. Therefore, using Mo.ipc to 
characterize the subhalo's mass allows us to consider subha- 
los that are more compact than standard virialized subhalos 
and makes it easy to apply our results to more exotic forms of 
dark matter substructure. 

To relate Qa to the virial mass of the subhalo, we have to 
specify the full density profile. If y = 2.0, we will assume 
that the subhalo is a truncated SIS so that Eq. ( fTSl ) holds out 
to the truncation radius of the subhalo. In this case, Eq. ( l23T l 
tells us that 6a - S^^, and we can use Eq. ( fTOl l to evaluate 6a- 
lfy + 2, we will assume that the subhalo's full density profile 
prior to any tidal stripping was 
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(25) 



which reduces to Eq. ( fTST i for r <s ro. In this case, the 
virial mass does not uniquely determine 6a, and we also 
have to specify the subhalo's concentration. We define the 
concentration as c = /?vii7'"-2, where r_2 is the radius at 
which d\npld\vir - -2. For the profile given by Eq. 
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r_2 = (2 - 7)ro. It follows that 
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(26) 

We see that ro is larger than 5 pc for the subhalos we consider 
(Mvir > 10"^ Mq and y > 1 .2), and we expect that local subha- 
los will be stripped to much smaller radii by encounters with 
stars. 

We can now derive how 6a depends on the subhalo's virial 
mass and concentration. Recall from Eq. ( l22l i that 0„ oc por^. 
For the profile given by Eq. ( l25T l, this factor is related to the 
subhalo's concentration and virial mass through 
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where B[z; a, b] is the incomplete Beta function. In Ap- 
pendix IB. 21 we use the findings of the Aquarius simulation 
JSpringel et al. 2008) to derive a relationship between the con- 
centration of local subhalos and their virial mass: 



c = 94 



Mv, 



10'' Mp 



-0.067 



(28) 



and we use this relation to determine the subhalo concentra- 
tion throughout this investigation. 

Figure U shows the paths taken by the star's image as the 
center of a subhalo passes 1 arcsec below the star's position 
for several values of y. In this figure, d^ - 50 pc, ds. - 5 
kpc, vt = 200 km/s, and the subhalo lens has a virial mass of 
5x10^ Mq. We see that the image path depends very strongly 
on y. If y ^ 1, the motion is nearly linear, and the image 
moves very little and very slowly. As -y increases, an arc ap- 
pears in the image path; there is now sufficient mass enclosed 
in the inner arcsecond to cause an observable vertical deflec- 
tion when the subhalo passes beneath the star. Increasing y 
also increases the motion of the image in a given time period, 
and the acceleration of the image as the subhalo approaches 
the point of closest approach becomes apparent. As the sub- 
halo center passes from the left to the right of star, the star's 
image will jump from right to left; since the image moves 
very slowly in the years before and after the crossing of the 
subhalo, this shift in the image's position off'ers the best hope 
for detection. 

3. OBSERVING STRATEGY 

The image motions shown in Fig. |4] suggest a simple detec- 
tion strategy for subhalo lensing, illustrated in Fig.|5] A typi- 
cal high-precision astrometric search program operates for up 
to 10 years. We start with an initial few-year calibration pe- 
riod, during which we assume that the star is relatively far 
from the lens. During this calibration period we 1) measure 
the star's intrinsic proper motion, parallax, and starting po- 
sition and 2) search for binary stars or other false positives. 
Stars that show significant acceleration during the calibration 
period probably have binary companions, and we reject them 
from the rest of the search. We follow the calibration period 
with a several-year detection run, when we essentially wait 
for a subhalo lens to come close to one of our target stars and 
induce significant lensing. 
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Fig. 5. — The subhalo lens detection scheme. The dashed line shows the 
trajectory of an image induced over 6 years by a subhalo lens with the same 
properties as in Fig.|4]and y = 1.8. The image motion during the calibration 
period and the detection run are labeled. The dotted hne shows the direction 
of the proper motion fitted during the calibration period, and the difference 
between the measured trajectory and the proper motion prediction is shown 
by the solid curve. For claiity, the intrinsic proper motion and parallax of the 
source are not shown. 

The timescale of the calibration period is important, as it 
needs to be long enough both to detect binary systems and to 
obtain a robust parallax and proper motion measurement. In- 
creasing the calibration period length improves the predicted- 
position accuracy, but it also takes time away from the de- 
tection run, reducing the probability of observing a lensing 
event. The calibration period must be at least one year long 
to obtain a secure parallax, and we suggest two years as a 
sensible length to ensure an accurate parallax and proper mo- 
tion measurement. The length of the detection run is then set 
by the duration of the high-precision astrometric campaign. 
In this analysis, we assume a four-year detection run, imply- 
ing a total of six years of observations, which is close to the 
expected mission lengths of Gaia and SIM and a reasonable 
length for long-term ground-based surveys. We leave a full 
discussion of the optimal observing strategy to future work, 
as the details of the observing scenario will depend on both 
the particular astrometric technique and the time available for 
the observations. 

For the purposes of this paper, we adopt a simple measure- 
ment of the astrometric signal from a lensing event: 



' ' cpoch.s 

. j] (Xim - XpiY + {Ynii - YpiY 



(29) 



where Xnij and Fm, give the 2D measured position of the 
star at each epoch, and Xp-, and Ypi give the 2D position of 
the star predicted from the proper motion, parallax and po- 
sition determine d during the c alibration period. This calcu- 
lation (following lGaudi & Bloo m 2005) essentially measures 
the total displacement of the star from its expected position 
over the course of the measurements. The signal-to-noise ra- 
tio (S/N) is simply S divided by the astrometric uncertainty 
per ID datapoint (cr); cr includes contributions from the in- 
strument's intrinsic uncertainty per datapoint as well as un- 
certainty in the star's proper motion, parallax and position. 
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The uncertainty in the star's predicted position grows in time 
due to the proper-motion uncertainty, and the parallax uncer- 
tainty's effect on the 2D measurements varies in direction and 
magnitude across the sky. 

This calculation gives a conservative estimate of the S/N of 
a possible lensing detection. The displacements induced by 
lensing are all in approximately the same direction, however, 
and the effective S/N would likely be improved by fitting a 
model to the data that accounts for these correlations. We 
leave such enhancements for future work. 

3.1. False positive removal 

In addition to simply detecting a subhalo lensing signal, we 
must also distinguish it from other astrometric signals. Sub- 
halo lensing signals take place over months, do not repeat, oc- 
cur without a visible lensing source, and have a near-unique 
trajectory. The relatively short event timescales ensure that 
our measurements are only sensitive matter structures on the 
spatial scales we consider here. The events' other properties 
can be used to remove false positives, such as those generated 
by stellar microlensing and motion in binary systems. 

Astrometric microlensing by a passing point-like lens 
moves the image centroid along an elliptical path that be- 
comes more circular as the impact parameter increases 
dWalker"1995; 'Paczv nskilfTogsl [19981: |MiraLda-Escude"1996'; 
[Gould 2000). The image path can become more complicated 
if the lens has a smal l finite extent and is opaque ( Ta kahashil 
l2003t[Lee et alJ[2010h . In all cases, however, the image com- 
pletes its orbit on observable timescales, quickly approaching 
its true position as the compact object moves further away. We 
saw in the previous section that subhalos produce a radically 
different lensing signature; after the passage of the subhalo 
center, the image moves very slowly and does not approach its 
true position until the edge of the subhalo passes by the star 
many decades later Lensing by subhalos is therefore readily 
distinguishable from lensing by compact objects. 

The astrometric signal of a stellar binary is easily distin- 
guished from subhalo lensing signals in most cases simply 
because the binary system produces a repeating signal. Al- 
most all long-period systems that do not produce repeating 
signals in our dataset will be removed by our requirement that 
the star does not accelerate during the calibration period. Fur- 
thermore, roughly circular binary systems induce a very dif- 
ferent astrometric signal from subhalo lenses. 

Rare highly eccentric systems with periods much longer 
than our observation length (or even very close unbound stel- 
lar encounters) can produce a short-timescale astrometric sig- 
nal during periastron, with little signal throughout much of 
rest of the orbit. The trajectory is similar to a lensing sig- 
nature (with an additional very large radial velocity signal). 
However, simple simulations of such systems show that no 
Keplerian orbit (of any eccentricity < 0.999) that produces a 
lensing-like signal can avoid producing detectable accelera- 
tion both during the calibration period and after the putative 
lensing event. 

3.2. Final subhalo lens confirmation 

The ultimate test of a candidate subhalo is a prediction of 
lensing. We will show in the next section that detectable dark 
matter subhalos are probably within 100 pc of our location; 
while it is possible to detect a more distant subhalo, the sub- 
halo would have to be massive (Mvir ^10^ Mq), and we ex- 
pect such objects to rarely pass between us and a target star. 



Detectable dark matter subhalos are therefore likely to have 
proper motions greater than half an arcsecond per year, and 
they will astrometrically affect all stars within several arcsec- 
onds of the subhalo center If a halo is detected, its path can 
be predicted (albeit initially at low precision), and a catalog 
prepared of stars that are likely to be affected by lensing over 
the next few years. Intensive astrometric monitoring of faint 
stars in the field could then provide a fairly rapid confirmation 
of the existence of the subhalo lens. 

4. CROSS SECTIONS FOR ASTROMETRIC LENSING BY SUBHALOS 
4.1. Signal Calculations 

As described in Section [3] we calculate the lensing signal 
by taking the square root of the sum of the squared differ- 
ences between the star's position on the sky and the position 
predicted by the proper motion and parallax measured during 
the calibration period [Eq. (|29]|1. At each epoch, the differ- 
ence between the measured image position and the predicted 
position is proportional to the deflection angle a. This linear- 
ity implies that we may use any vector fj that is proportional 
to a to calculate the signal; we just have to multiply the re- 
sulting signal by a/rj to obtain the physical signal that would 
be measured during a lensing event. We use this technique to 
separate the signal's dependence on the geometry of the lens- 
ing scenario from its dependence on the properties of the lens. 
As we describe in detail below, we can then calculate the geo- 
metrical signal once and then use that result to determine the 
signal for any lens. 

Following our convention, we work in a coordinate system 
where the subhalo's transverse velocity vj lies along the x- 
axis. In this case, /3y is the fixed impact parameter, and /3xfl 
specifies the initial position of the lens. It is useful to define 



If = and p = 



VT^obs 



t^'T^obs 



(30) 



where /3x,o and /?,. are in radians, and fobs is the length of the 
observation (not including the calibration period). If we de- 
fine A^ to be the angular distance, in radians, traversed by the 
lens during the observational period, we can easily interpret tp 
andjS. The normalized impact parameter's - /3y/A/3, while the 
phase (f = yS.v_()/Aj6 specifies the location of the point of clos- 
est approach to the star along the subhalo's path. Since we are 
only interested in cases where the subhalo center passes by the 
star during the observational period, we constrain < ip < I. 
With these definitions Eq. ( ITTT i may be rewritten as 



a ^Oai -^-^ I r]((p,^, t/tohs, y). 



(31) 



which allows us to separate the geometry of the lens-star sys- 
tem (i.e. ip and y6) from all of the lens characteristics apart 
from y. From Eq. (|2T] ). we see that 
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We use if to calculate the signal of a lensing event. This 
is advantageous because the resulting signal is independent 
of 6a, di, and vt; we call this signal Sg (for "geometrical sig- 
nal") because it depends only on if, fi, and y. The calculation 
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Fig. 6. — Left: The signal 5^ as a function of phase ((f) and impact parameter yS = fiy/A/} for a two-year calibration period and a four-year detection run. Note 
the decrease in signal toward larger phases, where the image only partially traverses its lensing path during the observational period. The decrease in 5^ at 
the smallest phases and impact pai'ameters corresponds to lenses that start to produce large image motion during the cahbration period, which is then partially 
subtracted out by the proper motion removal during the detection run. At larger values of ,§ the apparent motions are still relatively large (see Fig.|3] for instance), 
but these motions are difficult to distinguish from the star's proper motion, leading to the decrease in Sg. Right: The area A,, in the if - /3 plane that produces a 
geometrical signal that is larger than a given value for Sg for a variety of y values. Note that the normalized area goes to zero at a value of Sg that is dependent 
on y. For Sg values much smaller than this cut-off, Ag is proportional to (5^)'', where the index p depends on y. 

of Sg takes into account the subtraction of the proper motion infinity as the impact parameter's goes to zero. Consequently, 



and parallax measured during the calibration period, includ- 
ing any apparent motion generated by lensing, as shown in 
Fig-E] Figure |6] shows Sg{(p,^) for -y = 1.5 and y = 1.8. We 
see that Sg decreases with increasing impact parameter ;§ and 
increases with increasing y, which is not surprising given the 
image paths shown in Fig. |4] We also see a preference for 
geometries in which the subhalo passes by the star earlier in 
the observational period; the signal is enhanced because there 
are more epochs after the shift in the star's position when the 
subhalo center passes the star 
To relate the physical signal S to the geometrical signal Sg 
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This relation completes the procedure for determining if a star 
is detectably lensed by a given subhalo. Given a specific lens 
and a minimal detectable value for the signal S^in, Eq. ( |34] | 
may be inverted to obtain the corresponding minimal value 
for Sg. We then determine the area Ag of the region in the 
ip - P plane (with < i^ < 1 and p > 0) that produces a 
geometrical signal that exceeds this minimal value for 5^. Fi- 
nally, we convert Ag to a physical area on the sky that gives 
S > Smin, and we call this area the cross-section A(5'niin); in 
square radians we have 



A(5™i„) = 2|^| Ag(Sg). 



(35) 



The factor of 2 accounts for the fact that the stars both above 
(y§ > 0) and below (j6 < 0) the lens are equally deflected, and 
the two factors of {vjtobs/dL) follow from the definitions of tp 
andyg [see Eq. (poll]. 

The geometrical area functions Ag(Sg) are shown in the 
right panel of Fig. |6]for several value of y. As indicated by 
the Sg(sp,fi) surfaces shown in the left panel of Fig. |6] the area 
Ag that produces a geometrical signal larger than a specific 
Sg value depends strongly on y and decreases rapidly with in- 
creasing Sg. The left panel also shows that Sg does not go to 



there is a maximal value of Sg that is obtainable for each value 
of y. At this value of 5^, the area Ag goes to zero, as shown in 
the right panel of Fig. |6] This maximal obtainable value of Sg 
implies that there is a minimum subhalo mass that is capable 
of generating a detectable signal, as we wiU see in the next 
subsection. For values of Sg that are smaller than half of the 
maximal possible Sg value, Ag oc (SgY, as shown in the right 
panel of Fig.|6] For the four y values we tested, we found that 
p - -l/y. This simple power-law behavior will be shared by 
the physical lensing cross sections. 

4.2. Properties of the Lensing Cross Sections 

The basic shape of the lensing cross sections can be de- 
duced from the left panel of Fig. |6] This figure shows that 
Sg does not depend strongly on the phase ip when Sg is 
smaller than the maximum obtainable geometrical signal. For 
Smm values that correspond to these values of Sg, the lens- 
ing cross sections are rectangular The width of the rectangle 
is given by the change in the lens position during the obser- 
vational period (A/?) because we only consider lensing sce- 
narios with < (/? < 1. The length of the rectangle is two 
times the maximum impact parameter that a star may have 
and still be detectably lensed. The impact parameters that 
lie within the lensing cross section are therefore bounded 
by 1/3^1 < A(5'min)/(2A/S). For a lens at a distance of 50 pc 
with vt = 200 km/s, Eq. (|2]i tells us that A/3 ^ 3 arcsec- 
onds for a four-year observational period, so the maximum 
impact parameter that produces a lensing signal S > Smin is 
Py - MS min)/ (6 aicsec). 

Figure |2] shows how the lensing cross section depends on 
the mass of the subhalo for three values of ^min: 5 yuas, 20 yuas, 
and 50 /ias. We characterize the mass of the subhalo in two 
ways, as discussed in Section |23] In the top row, we show the 
cross section as a function of the mass enclosed in the inner 
0.1 pc of the subhalo (Mo.ipc). Recall that truncating the sub- 
halo density profile at a radius of 0. 1 pc does not significantly 
alter the lensing signal, which implies that Mo.ipc directly de- 
termines the lensing signature. In the bottom row, we show 
how the lensing cross section depends subhalo's virial mass. 
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Fig. 7. — The area around a subhalo that will produce an astrometric signal greater than S„,i„ as a function of subhalo mass, with d]^ = 50 pc, rfs = 5 kpc. and 
vt = 200 km/s. Top: The mass of the subhalo is defined as the mass contained within a radius of 0. 1 pc from its center. Truncating the subhalo density profile at 
a radius of 0. 1 pc does not aff'ect its lensing signature. Bottom: The mass of the subhalo is defined as 0. 1 % of its virial mass. The area curves for y = 1 .2 are not 
shown because the subhalo must have Myir > 10** Mq to generate a signal of at least 5 fias if y = 1.2. 



We assume that 99.9% of the subhalo's virial mass has been 
lost due to tidal stripping and take 0.00 IMyh as the present- 
day mass of the subhalo. 

From Fig. |7] we can see how the inner slope of the den- 
sity profile determines the strength of the lens signature; in 
all cases, the lensing cross section decreases sharply with de- 
creasing -y. The dependence on -y is stronger when the virial 
mass is used to define the mass of subhalo because a shal- 
lower density profile requires a larger virial mass to get the 
same mass within a given radius. The bottom row of Fig. Q 
indicates that subhalos with y < 1.5 are not capable of gener- 
ating easily detectible astrometric lensing signatures, and sub- 
halos with 7 = 1 .5 are only detectible if a star passes within 
a couple of arcseconds of the subhalo center The situation is 
far more promising for subhalos with y > 1 .8; in this case, an 
intermediate-mass subhalo could produce a signal of up to 50 
yuas if a star passes within 10 arcseconds of the subhalo cen- 
ter Moreover, we see that small subhalos (M, < 1000 Mq) are 
so concentrated that decreasing the inner slope of the density 
profile from y = 2 to y = 1.8 does not significantly change 
the astrometric lensing signal. Finally, the top row of Fig. Q 
shows that if the inner regions of the subhalos are denser than 
predicted by their virial masses, subhalos with shallower den- 
sity profiles are capable of producing detectible signals. 

For Smin values that correspond to geometrical signals well 
below the maximal possible value for Sg (the value of Sg at 
which Ag = 0), the cross section has a simple dependence 
on S„iia'- A is proportional to (5'mm)'', with p - -l/y for the 
four values of y we consider This simple power law is di- 



rectly inherited from the geometrical area functions Ag(Sg) 
discussed in the previous section and shown in the right panel 
of Fig. |6] Consequently, the power index p is independent 
of dL,ds,VT:, and subhalo mass, and it is even independent 
of whether Mo ipc or Myii is used to characterize the subhalo 
mass. As Smin increases toward the maximal possible value, 
the cross section decreases faster than (Smin)'' and rapidly goes 
to zero when the maximal possible value for 5min is reached. 

4.3. Lensing Event Rates 

The lensing cross sections computed in the previous sec- 
tions may be combined with a subhalo number density to yield 
a probability that any given star on the sky will be detectably 
lensed during the observational period. In this section, we 
will compute these probabilities using three candidate sub- 
halo number densities. The first two calculations will assume 
that all the subhalos have the same mass; we will assume 
that a fraction / of the local dark matter halo is composed 
of subhalos with mass M and radii of R - 0.1 pc, and then 
we will assume that the local dark matter halo was once in 
subhalos with virial mass Mvir. We take the local dark matter 
density pdm to be 0.4 GeV cm"-', which implies that there is 
3.5 X 10** Mq of dark matter within 2 kpc. Finally, we will use 
a subhalo mass function derived from the Aquarius simula- 
tions (Springel et al. 2008). Throughout this section, we will 
use vj - 200 km/s and fobs - 4 years. 

We compute the lensing probabilities by summing the indi- 
vidual cross sections for all subhalos with d]^ < ds for some 
fixed value of ds. We choose ds to be small enough that we 
may neglect the spatial variation in the subhalo number den- 
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Fig. 8. — The probability that a star is detectably lensed, with signal 5 > 5„,jn, if all subhalos have the same mass M within the inner 0. 1 pc. The local number 
density of the subhalos is /pjn, /M, where p^^^ is the local dark matter density, and / is the fraction of the dark matter contained within 0. 1 pc of a subhalo center. 
The lensing probability is proportional to /, and / = 1 is shown here to illustrate the maximal possible lensing probabihty. 



sity within this sphere, and we assume that the subhalos are 
isotropically distributed. From Eq. (l35T l, we see that the lens- 



ing cross section for an individual halo is A oc d^Ag, and A 



depends on d]^ only through the (1 - dL/ds) factor in 9a. 
follows that the total cross section is 



It 



Itot = ^Jr (vifohs) ds nsub(M) 



Jo 



9a(x = dL/ds) {v-rt, 



1-0 



2-y 



dx, (36) 



where nsub(A/) is the number density of subhalos with mass 
M. Although (is will typically vary from less than 1 kpc to 
5 kpc for stars in a high-precision astrometric survey, we sim- 
plify the calculation by assuming that all the monitored stars 
are 2 kpc from us. Since Atot increases linearly with ds, taking 
ds. - 2 kpc effectively averages over stars that are uniformly 
distributed over 1 kpc < t/s < 3 kpc. 

We will see that the resulting total cross section Atot is much 
smaller than the total area of the sky As^y- It is therefore 
highly unlikely that the cross section for lensing by one sub- 
halo will overlap with the cross section associated with a dif- 
ferent subhalo, and we may consider Atot to be the total area 
on the sky in which a star would be detectably lensed. Fur- 
thermore, since the individual subhalo lensing cross sections 
are less than 0.1 square arcminute (see Fig.|7]i, the probability 
of having multiple stars within the cross section of a particu- 
lar subhalo is low enough that we may consider each star to 
be an independent sample of the sky. In this case, we may in- 
terpret the fraction Atot/Asky as the probability that any single 
star is detectably lensed by a subhalo. A subhalo in the Galac- 
tic plane is more likely to be detected than a subhalo near the 
Galactic pole due to the higher density of target stars in the 
plane, but, since the subhalos are isotropically distributed and 
can only lens one star apiece, the average lensing probability 
(Atot/Asky) is individually applicable to each star in the sky. 

4.3.1. Event Rates: Mono-mass subhalos 

We first consider cases where all the subhalos have the same 
mass Mo.ipc within a radius of 0. 1 pc. The local subhalo num- 
ber density is then Wsub - fPdm/Mo.ipc, where / is the fraction 
of the dark matter that is contained in the inner 0. 1 pc of the 
subhalos; / < 1 could mean that there is a smooth component 
of the local dark matter distribution, or it could mean that the 
subhalos' radii are larger than 0.1 pc and their actual masses 
are Mo.ipc//- In Fig- El we show the lensing probability with 
/ = 1 for three values of ^min: 5 //as, 20 /zas, and 50 /zas. We 



see that the lensing probability is highest for each value of 
y if the subhalo mass is just slightly larger than the minimum 
mass required to generate a sufficiently large signal- Although 
more massive subhalos have larger cross sections, as shown in 
Fig- 13 the dependence of A on Mq ipc is not steep enough to 
compensate for the diminishing number density of subhalos 
as Mo.ipc increases- 

If we instead assume that all of the local dark matter was 
originally in subhalos with virial mass Mvii and that the cen- 
tral regions of these subhalos survive to the present day, then 
the subhalo number density is Wsub - Pdm/Mvir- As in the pre- 
vious case, the (l/Mvii) factor in «sub implies that the lensing 
probability will peak near smallest value of Mvir that is ca- 
pable of generating a signal- For S^in - 5yuas, the lensing 
probability peaks at Mvir - IO'^Mq for 7 = 2 and y = 1-8, 
and Mvii - 10*^ Mq for y - 1.5, as predicted by Fig.|7] Since 
the cross section for these smaller subhalos does not change 
much between y = 1.8 and y - 2.0, the lensing probabilities 
for these two cases are very similar and they peak at probabil- 
ities of 5 X 10"'" and 7 x 10"'", respectively. The y = 1.5 case 
peaks a far lower probability of 1 x 10"'^. It is not surprising 
that these numbers are about four orders of magnitude lower 
than the peak probabilities in Fig. [HI from Fig.|7] we see that 
A(Mo.ipc) ~ A(Mvir) for Mo.ipc ~ lO^^'Mvir and y = 2, so by 
setting «sub - Pdm/Mvii-, we are effectively setting / ~ lO"'' in 
Fig- El 

When we assume that all the subhalos have the same mass, 
then Atot has the same simple dependence on Smin as A; for 
Sg values that are small compared to the maximum accessible 
value, Atot i s proportional to {SntinV, with p = -1/y as in 
Section 14-21 If we instead consider a subhalo mass function 
and integrate over subhalo masses, the dependence of Atot on 
Smin changes because there is no longer a single value of Smin 
that corresponds to the value of Sg at which Ag = 0. 

4-3-2- Event Rates: subhalo mass function 

To evaluate the lensing probability with a range of subhalo 
masses, we will u se the local subhalo mass function derived 
in AppendixIB-llfrom the results of the Aquarius simulations 
JSpringel et al.ll2008l) : 



dwc 



dMv 



= 2x10" 



Mvir 

Mo 



-1.9 



(37) 



As described in Appendix lB.il this subhalo mass function is 
applicable within a few kpc of the Sun, and we have assumed 
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that the subhalos in this region lose 99% of their virial mass 
due to tidal stripping by the smooth component of the dark 
matter halo and other subhalos (stars are not included). We 
evaluate Atot by replacing Wsub in Eq. ( l36b with dnsub/dMyir 
and integrating over Myir from M^in to M^^^. For each value 
of y and S^m, there is a minimal value of Mvir needed to make 
Ag nonzero; this minimal virial mass is always larger than 
10 M0, so we set M^in - 10 Mq. We choose M^ax such that 
the expectation value for the number of subhalos with Myii > 
Mmax within 2 kpc is greater than one; from Eq. (|37] |. we have 
A^max = 3 X 10^ Mq. Extending the integral to larger values 
of Mvir changes the meaning of Atot; it is no longer a sum 
of cross sections for the subhalos expected to be within 2 kpc. 
Instead, Atot would also include the cross sections for subhalos 
that we do not expect to find within 2 kpc, multiplied by the 
probability that such a halo is in this volume. 

The value of M^^^ determines the shape of the AtotiSmin) 
function because it determines the value of S^in at which Atot 
goes to zero. If Sg{Smm,Mm^,,) is greater than the maximum 
reachable value for Sg, then no subhalo with Mvir < ^max can 
generate a signal and Atot - 0. For -y - 1.2, 1.5, 1.8, and 2.0, 
this maximum value for Smm is 0.46, 10, 250, and 750 yuas, 
respectively. If Smin is much less than these upper limits, then 
Atot(5n,in) is roughly a power law, and the lensing probability 
is approximately 



Atot 

A sky 

Atot 
-A sky 



7=1.8 



: 8.7x10-12^1!!^ 
\5/ias 

= 1.3xlO-"(^^ 
y=2.o \5^as 



for^min < SOyuas, (38) 



for^tnin < 200//as.(39) 



For larger values of 5min, Atot(5tnin) decreases faster than these 
expressions and quickly goes to zero at the values listed 
above. Ify = 1.5, then Atot/Asky = 1.4xl0-'2for5niin - lyuas, 
and there is no power-law behavior between S^in - 1 yuas and 
the zero point at 5„,in = 10/ias. 

Including the probabilities that larger subhalos are present 
within 2 kpc of our location does not significantly affect the 
total lensing cross section. If we set M^ax = lO'^M©, then 
the dependence of Atot on S min is slightly shallower than the 
power laws given in Eqs. (l38T l and (|39] |. but the differences 
are not large. For S'lni,, < 20;uas, including subhalos with 
3 X 10^ M0 < Mvii- < 10'" Mq increases the lensing probabil- 
ity by less than 20%. At the largest Sj„i„ values described by 
Eqs. (l38T l and ( |39] |, extending the integral to these larger sub- 
halo masses increases Atot by a factor of three. The inclusion 
of larger subhalos has a more pronounced impact on the value 
of Atot for Smin > 0.46, 10, 250, and 750 yuas for y = 1-2, 1.5, 
1.8, and 2.0, respectively. The larger value of M^ax implies 
that the total cross section does not go to zero at these values 
of 5n,in as it did when we set Mmax = 3 x 10^ Mq. This is not 
an important change, however, because the lensing probabil- 
ity is less than 10^ for these large values of 5 min. Moreover, 
even if it were possible to monitor far more than 10'^ target 
stars, lensing events would only be observed if a subhalo with 
Mvir > 3 X 10^ Mq lies between us and the target stars. 

5. DETECTION PROSPECTS 

As discussed in the introduction, it is only recently that as- 
trometric capabilities have advanced to levels where the de- 
tection of subhalo lenses becomes possible. In this section 
we use the results derived above to evaluate current and fu- 
ture astrometric subhalo search techniques. We consider two 



scenarios: 1) a large-area search for subhalo lenses and 2) a 
confirmation of a subhalo suspected on the basis of other de- 
tection methods. 

5.1. Achievable Smin 

We start by evaluating the achievable S^in for an astromet- 
ric survey. For a particular astrometric observing strategy, 
both the number of epochs and the astrometric precision affect 
Smin- For the purposes of this discussion we calculate the best 
Smin for each technique using the six-year observation setup 
described in Section [3] We include the single epoch instru- 
mental uncertainty, as well as a detailed Monte Carlo simu- 
lation of the extra per-epoch uncertainties introduced by the 
measurement and subtraction of the target star's calibration- 
period position, proper motion and parallax. We marginalize 
over a wide range of possible parallaxes and proper motions, 
as well as the full range of possible sky positions. 

For a particular observational setup, the simulations pro- 
duce a scaling factor between the instrumental astrometric 
uncertainty per epoch (cri„st) and the final total uncertainty 
per datapoint (cr) used to calculate S/N. For the observational 
setup described in Section [3] cr = 1.47crinst. For a typical 
S/N - 3 detection, an instrument's S^in is then 4.4crinst. We 
note that the summed nature of Smin means that astrometric 
displacements smaller than crinst are indeed detectable in this 
scheme. 

5.2. Large-Area Surveys 

A practical all-sky search for subhalo lenses requires 
enough stars that the lensing probabilities described above 
lead to a significant probability of detection. For example, 
following the left panel of Fig. [8] if all subhalo lenses are sin- 
gular isothermal spheres and have a Mo.ipc of 2 Mq, and our 
survey is capable of detecting 5',v,in=5 fias, we need to survey 
~5/-' X 10^ stars to have a good chance of detecting a sub- 
halo lens, where / is the fraction of local dark matter that is 
contained within 0. 1 pc of a subhalo center Alternately, at an 
Smin of 50 jt/as we need to survey ~3/-' x 10* stars. 

Ground-based all-sky surveys ar e currently limited to milli- 
arcsecond precisions at best (e.g.. Ilvezic et alJl2008l) . so we 
do not consider them further here. The Gaia mission will 
have an astrometric uncertainty of crinst - 35 microarcsecond 
per epoch at my - 12 for its best targets, with an average of 
83 epochs per targefl- From a search of the USNO-Bl cata- 
log, we find that ~ 5 x 10* stars are covered at this precision 
(lMonetetalJl2003h . For these tai-gets. Gala's 5,™nis -260 
//as because coverage of 10* stars requires S/N~5 to avoid 
false positives; this precision is too low to detect intermediate- 
mass subhalos. Orders of magnitude more stars are cov- 
ered by Gaia at lower precision, but an extremely large (and 
very unlikely) subhalo mass would be required to produce de- 
tectable lensing at those precisions. Although there remains 
a small probability that Gaia could see a lens, we conclude 
that Gala's astrometric precision is probably insufficient for 
a useful all-sky subhalo search. That said, it is prudent to at 
least attempt a subhalo lens search using the Gaia dataset, as 
the data will be available and can be readily searched for such 
signals. 

For the N-bod y simu lation-based lensing probabilities dis- 
cussed in section lT.3.21 where typical lensing probabilities are 
~ 10"" at Smin = 5 fj.as , much improved instrumentation ca- 
pabilities would be required for detection. The only currently 

'' [http : //www ■ r ssd ■ esa ■ int/index . php?pro j ect $=$GAIA&page$=$Science_Per 



Local DM Subhalo Lensing 



13 



planned instrument capable of reaching S'niin = 5 //as routinely 
is SIM, but SIM will at most cover tens of thousands of tar- 
gets during its lifetime d Unwin et alj 120081) . If the subhalo 
mass functions derived from simulations are correct, a SIM- 
precision all-sky search would have to cover 10" targets to 
have a good chance of making a detection. This capability 
will most certainly not be available in the near future. How- 
ever, it is worth stating that a SIM-precision mission cover- 
ing 10** targets (a possible capability for a next-generation all- 
sky astrometric survey), would place unique constraints on 
the subhalo mass function. In particular, such a survey could 
usefully evaluate if the simulations under-predict the subhalo 
mass function or if the subhalos are more dense than expected 
(the scenario in Fig. |8]i. 

Similar conclusions can be drawn for blind searches of sub- 
areas of the sky: current and planned astrometric capabili- 
ties are insufficient for a large-area survey for subhalo lenses. 
However, if a local subhalo is suspected on the basis of other 
detection methods, targeted surveys are capable of either de- 
tecting the lens directly or stringently constraining its proper- 
ties. 

5.3. Targeted Observations 

If we have some idea of where a lens might be, searching 
for that lens becomes much easier For examp le, it has been 
suggested that the Fermi Large Area Telescope dAtwood et al.l 
[2009) could detect subhalos in gamma-ray emission (e.g., 
lUegal-Gaskins 2008; Ando 2009; Bucklev & Hooijer 2010). 
Fermi's first point source catalog (Abdoetal. 2010) con- 
tains a large number of unidentified sources that could be 
subhalos capable of pro ducing detectable lensing signals 
jBucklev & Hoopeil20I0l) . Sources in the Fermi point-source 
catalog are localized to 6 arcminutes (median; 95% confi- 
dence) or even 1.5 arcminutes (best 50 sources; 95% con- 
fidence). The lensing search therefore requires coverage of 
only 0.01 sq. degrees or less. Furthermore, since Fermi has 
many plausible sources, we can pick the targets with the best 
likelihood of detection, such as sources close to the Galactic 
plane with many astrometric target stars. The aim here is to 
place the best possible limits on the lens properties (with the 
possibility of an actual detection), and even current techniques 
(reaching S^m - 50 //as ) could place useful limits on the lens 
properties. If astrometric lensing is detected around a gamma- 
ray source, then the magnitude of the deflection provides a 
measurement of the subhalo's central density; the shape of 
the image's path provides information about the inner den- 
sity profile; and the rate of change in the image's position 
provides information about vj/c/l- If no lensing is detected, 
constraints could be placed on a combination of the object's 
distance, mass, and density profile. 

Since it is unlikely that we will know the exact position of 
a suspected subhalo, moderately wide-field astrometric tech- 
niques are favored for this type of search. From space, the 
Hubble Space Telescope (HST) has demonstrated ~1 mas pre- 
cision crowded-field astrometry in the cores of globular clus- 
ters using the Advanced Camera for Surve ys (ACS) instru- 
ment (e.g.. Anderson & van der Marell2010l) . However, much 
better precision has been demonstrated on arcminute scales 
from the ground. 

AO-equipped 5-lOn i telescopes routinely achieve 100 p.a s 
precision astrometry (ICameron et al.l l2b09t iLu et al.l 120091) . 
These techniques minimize systematics by observing in nar- 
row NIR bands, at consistent air masses, and with careful at- 
tention to other sources of systematic error. Such surveys re- 



quire only a field with several stars within an arcminute field 
of view (such as is common in the galactic plane) and a few 
minutes of observing time per field. Using these techniques, 
relatively small 2m-class telescopes equipped with low-cost 
adaptive optics systems can reach 50-100 //as precision in 
tens of r ninutes and can per f orm large, intensive astrometric 
surve ys (iBritton et all 120081: iBaranec et all 120091: iLaw et al.1 
120091) . The precision can be further improved; in the ab- 
sence of systematics, a lOm-class telescope is predicted to 
reach 10 //as in similar integration times (although currently 
systematics limit precisions to the ~100 //as level, develop- 
ment is continuing). On the same basis, a 30m-class telescope 
coul d reach few-/ias p recision in a few minutes over a small 
field (ICameron et al.l2 009). although the systematic errors are 
again likely to dominate such observations until technique im- 
provements are made. 

Current ground-based techniques with precisions of 
50-100 //as can detect ^mindown to ~200 //as with S/N=3. 
This precision may be enough to detect nearby large subha- 
los (Mvir > 5 X 10^ Ms or Mo.ipc > 400 Mq for y > 1.8) in 
the galactic plane with current instruments. If the astromet- 
ric accuracy is improved, current lOm-class telescopes could 
potentially reach Sm\n <50 //as in 10-minute observations. In 
crowded regions, such a system performing a long-term astro- 
metric survey could detect subhalos down to stripped masses 
(niM = 0.001) of -1000 Mq at -50 pc distances, while 30m- 
class telescopes could detect subhalos at least an order of 
magnitude smaller. 

SIM offers another possible route for subhalo lens confir- 
mation. The instrument can reach a best precision of 1 //as 
(and so a detectable ^min of ~5 //as ). Crowding limits require 
that SIM's target stars are separated by at least ~5 arcseconds 
from each other, and so the lensing area must subtend at least 
-25 sq. arcseconds. With these capabilities, SIM would be 
capable of confirming a lens 50 pc away down to a stripped 
virial mass of -100 Mq (for y =2.0 or 1.8), or equivalently 
MoApc of 10 Mq (Fig.|7]i. Because SIM is a pointed mission, 
it can target faint (mv=20) stars at the cost of observing time, 
making it relatively easy to obtain a sufficient sample of stars 
near to a putative subhalo position. Although it may not be 
possible to obtain 1 //as precision observations of all stars near 
to a suspected lens because of observation time constraints, 
we estimate that ^min of 4-18 //as would be obtainable in mod- 
est amounts (weeks) of SIM observation time. SIM observa- 
tions could thus confirm suspected subhalos near the galactic 
plane down to stripped subhalo masses (mbd - 0.001) of hun- 
dreds of solar masses or Mq ipc of tens of solar masses. 

6. SUMMARY AND CONCLUSIONS 

When a dark matter subhalo moves relative to a more dis- 
tant star, the star's apparent position changes due to gravi- 
tational lensing. By studying the image motion generated 
by subhalos with isothermal and NFW density profiles, we 
have determined that the change in the image's position is de- 
tectable only if the subhalo's center passes by the star with a 
small impact parameter (< 0.01 pc in the lens plane). There- 
fore, only the inner 0. 1 pc of a subhalo is relevant for astro- 
metric lensing, and we adopt a general power-law density pro- 
file in this region. We used the findings of the Aquarius sim- 
ulation (Springel et al. 2008) to derive a relationship between 
the concentration of local subhalos and their virial masses, 
which allows us to convert between a subhalo's virial mass 
and the mass enclosed within 0. 1 pc of the subhalo center. We 
found that the image paths due to lensing by a subhalo with 
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p cc r ~^ depend strongly on y, with cuspier profiles producing 
much larger deflections than shallower profiles. For y > 1.5, 
an intermediate-mass subhalo (Myir ^ 10'' Mq) within a kpc of 
the Sun can produce astrometric deflections that are detectable 
by current and near- future instruments. 

We have designed an observing strategy that can be used to 
detect subhalo lensing in data from typical astrometric sur- 
veys (such as Gaia, SIM, or ground-based methods). The 
setup makes use of the typical subhalo lensing signature: the 
image starts out essentially fixed in position, and as the sub- 
halo center passes by, the star rapidly moves to a new po- 
sition. The star's image then remains nearly stationary for 
the next several years. Under our observing scheme, a star's 
position, proper motion, and parallax is measured during an 
initial calibration period, and then the star is monitored over 
the next several years. We define the lensing signal as the dif- 
ference between the lensed trajectory and the path predicted 
by the measured proper motion and parallax. This strategy 
is immune to the most important false-positive possibilities: 
eccentric binary stars and point-source microlensing. 

The magnitude of the resulting astrometric signal depends 
on the impact parameter between the subhalo and the star For 
a given minimal signal required for detection, a given subhalo 
will detectably lens all stars within a certain area on the sky. 
We computed this cross section for lensing as a function of 
both the subhalo virial mass and the mass enclosed within 
0. 1 pc for several values of y. Combining these cross sec- 
tions with a local subhalo number density allows us to cal- 
culate the probability that a given star's image will be de- 
tectably deflected by a subhalo within a given observation 
time. To evaluate the subhalo lensing probability predicted 
by N-body simulations, we derive a mass function for local 
subhalos based on t he findings of the Aquarius simulation 
dSpringel et al.ll2008l) . 

Finally, we use these cross sections and event rates to eval- 
uate the detectability of subhalo lensing signatures. We con- 
sider two scenarios: 1) a large-area survey for subhalo lenses 



and 2) a confirmation of a subhalo suspected on the basis 
of other detection methods. We find that Gaia all-sky as- 
trometric measurements are close to being able to constrain 
subhalos with abnormally high central densities, as could 
arise if substructure formed very early ( Ricotti & Gould 20091 
Berezinsky et al. 2010) . A subhalo's astrometric lensing sig- 
nature would also be enhanced if it conta ins a black hole 
jBertone et all 120051: iRicotti & Gouldll2009h : the black hole 
would steepen the inner density profile and would add a point 
mass to the subhalo center, resulting in a distinctive astromet- 
ric lensing signature. Given these possibilities, it is certainly 
prudent to attempt a subhalo lensing search using Gaia, and 
we leave a thorough investigation of these more exotic sce- 
narios for future work. Unfortunately, a full-sky survey with 
much higher astrometric precision than Gaia is required to 
usefully constrain the dark matter subhalo mass function cur- 
rently predicted by N-body simulations. 

A targeted search for astrometric lensing by subhalos is far 
more promising; if a subhalo's presence is suspected by other 
means (for example, as a Fermi gamma-ray source) current 
and near-future ground based astrometry surveys are capa- 
ble of directly searching for the subhalo's lensing signature, 
down to stripped masses (O.OOlMyir) of ~1000 Mq at ~50 
pc distances. The SIM astrometric satellite could confirm 
suspected subhalos near the galactic plane even if the sub- 
halo is 1-2 orders of magnitude less massive or more distant. 
Fermi has already observed gamma-ray sources of unknown 
origin, and the possibility that the se sources are dark mat- 
ter subhalos has been investigated (iBuckley & Hoo per 2010). 
If Fermi detects a gamma-ray source that is consistent with 
a nearby intermediate-mass subhalo, then high-precision as- 
trometry could at a minimum constrain the object's distance, 
mass, and density profile, and it may even provide definitive 
confirmation for the detection of dark matter substructure. 
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APPENDIX 
STRONG LENSING AND TRUNCATION EFFECTS 

The lensing cross sections presented in Section W2\ confirm that we are firmly in the weak-lensing regime (a <k /? in Fig.[T]i. 
It follows from Eq. (l23T l that the condition a <K y6 is equivalent to the condition 9e <si j3. In our coordinate system, with the lens 
moving along the .^c-axis, the minimal value foryS is the impact parameter/?,,; we are in the weak-lensing regime only if 6e <k fSy 
We should excise the area with 6e > /3y from our cross sections because our solution to the lens equation is not valid in this 
region. If A/3 is the angular distance traversed by the lens during the observation period, then the area that should be excised is 
A V = 2A/3 X 6e- This area is much smaller than the total cross section for lensing, A = 2(Ap)^Ag from Eq. (l35T l. if A^ » Oe/A^. 
For all the subhalos that we consider (Mvii- < 10^ Mq, y < 2), S^/Ag < lOyuas. With this bound and vt > 5 km/s, A/3 » OE/Ag 
is satisfied if d]^ <k 1000 kpc. We therefore conclude that A , is an insignificant contribution to the lensing cross section and we 
need not exclude it. 

The cross section for lensing also tells us how far a star may be from the center of the subhalo and still be detectably deflected. 
As discussed in Section |231 we do not truncate the density profile when calculating the surface density of th e sub halo, and this 
approximation is valid only if ^ is much smaller than the truncation radius (R,) of the subhalo. In Section |23] we stated that 
this condition is safely satisfied if R, > 0.1 pc, and we will now verify that claim. We assume that we are only interested in 
subhalos with virial masses less than 10^ Mq and signals greater than 1 yuas. These restrictions define a lower bound on detectable 
values of the geometric signal; for y = {2.0, 1.8, 1.5, 1.2), we have Sg > {0.0020,0.0056,0.098, 1.22). At these small values, 
Sg is nearly independent of the phase ip, as illustrated in the left panel of Fig. |6] Since tp is confined to be between and 1, 
the maximal value of y6 that keeps Sg above these lower bounds is just the area Ag evaluated at the minimal value of Sg: for 
y = {2.0, 1.8, 1.5, 1.2), we have p < {37,28, 10,0.2). From the definition of yS, we see that ^y = vj^obsA ^"d so the maximal 

distance (in the lens plane) between the star and the subhalo center is ^^ax - vxfobs Vl+J^max- With fobs = 4 years, and vj - 200 
km/s, we have ^^^x = {0.03, 0.02, 0.008,0.0008) pc for y = {2.0, 1.8, 1.5, 1.2). If we take R,>0.1 pc, then ignoring the truncation 
of the density profile overestimates M2D(^max) by less than 20% for y > 1.5 and less than 40% for y = 1.2. Moreover, these are 
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very conservative estimates; for Mvir <k 10^ Mq or ^min s> 1 //as, the minimal values of Sg will be much larger, leading to smaller 
values of f max and less disparity between the truncated and infinite values of M2D(^max)- 

SUBHALO PROPERTIES 

We are interested in intermediate-mass subhalos (10 Mq < M, < 10* M©) that are located within a few kpc of the Sun. 
Fortunately, numerical simulations have recently reached the resolutions required to probe substructure within a few kpc of the 
Sun (ISpringel et al.ll2008l : iDiemand et al.ll2008l) . with minimum resolvable subhalo masses of ~ 10^ Mq. In this appendix, we 
will use the Aquarius simulation results presented by ISpringel et alj ( 120081) to derive a mass function and a concentration-mass 
relation for these local subhalos. 

The Aquarius simulation suite includes simulations of six galaxy-size dark matter halos with ~ 2 x 10*^ particles in each halo 
and a higher resolution simulation of one of these six halos (Aq-A) with 1 .4 x 10"^ particles. In the Aq-A simulation, each particle 
has a mass of 1712 Mq, making it possible to identify subhalos with masses greater than 10^ Mq. ISpringel et al.l (12008b defines 
the interior of the host halo as a sphere with a mean density that is 50 times the critical density; the radius of this sphere called 
rso, and the mass enclosed is Msq. The Aq-A halo has M^q - 2.5 x 10'^ M© and rso = 433.5 kpc. The mass of a subhalo (M,) 
is determined by the SUBFIND algorithm ( S pringel et alj|2001l) . which counts the number of gravitationally bound particles and 
then multiplies by the mass per particle to obtain M,. We will continue to use M, = WbdMvii to parameterize the effects of tidal 
stripping on the subhalo's mass. 

Subhalo Mass Function 
The subhalo mass function measured in the Aq-A halo for all subhalos with r < r^o is 

dA^ (MrY^'^ 

dM, \mol 

with fl() = 8.21 X IO'M'^q and mq - lO^^M^o- ISpringel et al.l (120081) also report that the subhalo number density has the same 
spatial dependence for all subhalo masses 10^ Mq < M, < 10'" M©: 



n{M, r) - no(M) exp 



2.. . . _^ 



a I\0.46r50 



(B2) 



with a - 0.678. To determine the function dno/dM,, we integrate Eq. (IB21 i over r < r^o and match the result to Eq. (IB lb . The 
resulting mass function is 

jtL.^5L_w''„p[_i|(^_y'_a|. ,B3) 

dM, 1.985 4, \ mo/ I aU0.46r5o/ J J 

The Aq-A halo is larger than the Milky Way's halo (e.g., iDehnen et al]l2006t iLi & Whitdr2008l; IXue et alJlIOOl iReid et alJ 
2009h. so we mu st use appropriate values of rso and Mso when evaluating Eq. (IB3l l. We use the density profile presented by 



Xue et al.l (120081) to derive approximate values of rso and Mgo for the Milky Way: Mgo - 9.5 x lO" Mq and rgo - 310 kpc. With 
these parameters, 

1 n _ 

exp 



dn 2.5 X 10-^ / M, ^ ''' 



dM, pc^^ Mq \ Ml 



iQ 



a ^ 140 kpc 



(B4) 



For r < 20 kpc, the subhalo number density is no longer strongly dependent on r, and it changes by only 7% as you move 2 kpc 
away from the solar radius (/?o - 8 kpc). We will neglect these small variations so that we may treat the local subhalo number 
density as isotropic. If we evaluate Eq. (IB4I) at r = 8 kpc, we find 

d« 3 X 10-^ / M, \"'" 



dM, pc^ Mq \ Mq 

We now need to convert this mass function for subhalo mass M, to a mass function for virial mass Myir. We will assume that 
all subhalos within 2 kpc of the Sun lose the same fraction of their mass due to tidal stripping. If M, = m^Myii, where m^ is 
constant, then we have 

dn 3xlO-^m,r/M™V''' ^gg^ 



dMvir pc^ Mq \ Mq 

There is great uncertainty surrounding the local value for mbd- IDiemand et al.l (l2007bh monitored subhalo mass evolution in the 
Via Lactea simulation, and unsurprisingly found that the fraction of mass lost due to tidal stripping increases closer to the center 
of the host halo. They found that subhalos in the region containing the inner sixth of the host halo mass lose roughly 80% of 
their mass between a redshift z ~ 2 and the present day. This sample contains subhalos t hat are far further from the host's center 
than the Sun, so we may consider mbd - 0.2 to be a rough upper bound. Meanwhile, van den Bosch et aLJ (1200 51) developed a 
semi-analytical model for tidal stripping and concluded that 0.001 < mb d $0.1 for all subhalos in a Milky- Way sized host, with 
most subhalos losing 99% of their original virial mass. Finally, Eq. ( |B6! l implies that the total mass in a sphere with radius 2 kpc 
that was once part of a subhalo with Myi, < 10^ Mq is 6m^''^'^ x 10^ Mq. This mass must be less than all the dark matter contained 
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in this sphere (3.5 x 10** Mq), so njbd ^ 0.001 on averag e. We adopt a middle-of-the-road value of m^ - 0.01 for local subhalos 
when evaluating the lensing event rates in Section !?. 3. 21 This value for ;«bd does not include stripping by stellar encounters; here 
M, is the subhalo mass measured in N-body simulations that do not include stars. 

Subhalo Concentrations 

Many methods for assignin g concentrati ons to dark matter halos have been proposed (e.g iBullock et al.l 120011; iNeto et aT] 
[2007; Duff y et 30120 08; Macc io et al.ll2008l) . but these models focus on isolated halos that are far more massive than the sub- 
halos we consider. Furthermore, numerical simulations indicate that subhalos nearer to the center of the host halo have higher 
concentrations t han both isolated h alos (Ghi gnaet al. 1998; Bullock et al. 2001) and subhalos in the outskirts of the host halo 
jDiemand et al.t2 007b; S pringel et a l. 2008; Diemand e t al. 2008). In light of this dis tinction, we adopt a c(Mvii) relation for local 
subhalos that is based on the findings of the Aquarius simulation (.Springel et al.ll2068,) . Since the virial radius of a subhalo is not 
easily measured, a different concentration parameter is often used to characterize the concentration of subhalos in simulations: 

2y^ 

6v ^ ttt"^^ (B7) 



max 



where Vmax is the maximum circular velocity within the subhalo, and r^ax is the distance from the subhalo center at which the 
circular velocity is maximized. The Aquarius team found that 6v depends on subhalo mass M, and distance from the host halo 
center r; when they average over all subhalos with M, > 3 x 10^ Mq, they find 



\-0.63 
^6/ 



6v = 3.8 X 10" , (B8) 

\kpc; 



and when they average over all subhalos, they fincQ 

= s 8 V 1 n4 1 

, 10^ M< 

with considerable scatter in both cases JSpringel et al.ll2068h . Inspired by these relations, we adopt a model 



5-8 X 10^ I ^^J , (B9) 



iZ-V'''i^L "°" 



and we use the position-dependent subhalo mass function derived in the previous subsection to compare this model with Eqs. 
( IB8l l and ( |B9t . Matching Eq. ( IB8l l gives A^^ = 2.4 x 10^, while matching Eq. ( |B9l ) gives A^^ = 1.5 x 10^. Since increasing the 
subhalo concentration enhances the lensing signal, we adopt the latter value to be conservative. 

Given a full (pre-stripped) density profile for the s ubhalo, it is possible to relate (5i/ to c = Rvulr-2, where r_2 is the radius at 
which c/lnp/c/lnr = -2 (e.g.. JDiemand et al]l2007bh . For the density profile given by Eq. ( |25] l. 



(5i/ = 



Pvii(Zv) 



Pcrit,0 



' g[ymax(y - 2); 3 -y,y- 2] 
Jmax/ B[c{y-2);3-y,y-l\ 



(BID 



where z,. is the redshift at which /?vii- is evaluated, pa-it.o is the present-day critical density, y^ax = '"max/'"-2, and B[z; a, b] is the 
incomplete Beta function. For 1 < y < 2, y^^^ ^2.1, and the function dyic) is not strongly dependent on y. Since we are 
interested in subhalos with 10"* M© < Mvh < 10^ Mq and r ^ 8 kpc, we only need to consider the range 10^ < 6v < 10^. In this 
range, 5v(c) is well-approximated by a simple power law: 



0.049 



Pvir(Zi.) 



c^"'. (B12) 



Pci-it.O 

We obtain a final expression for c(Mvii) by inverting Eq. (IB12I) and inserting Eq. ( IBIOI ). With Ng - 1.5 x 10^, we find 

/ _ , . x-0.37 / ,, X -0.067 / \-0.24 

, = 94,«,r^plil^ J^ U_] , (B13) 

>"* \4.6Mepc-V \106Mo/ \8kpc/ ^ ^ 

where we have defined mbd = M,/Mvii. Since c(Mvir) changes little for 0.01 < mbd < 1, we will take mbd - 1 when evaluating the 
concentration. The dependence of c(Mvir) on the subhalo's location is also fairly weak; for a fixed value of Myir, c decreases by 
only 10% between 8 kpc < r <13 kpc and increases by only 25% between 3 kpc < r <8 kpc. We conservatively take r = 8 kpc 
when evaluating the concentration of local subhalos. Finally, we not e that subhalos in simulat ions that include baryons tend to 
be more concentrated than subhalos in simulations without baryons dRomano-Diaz et ani2010l) . so it is possible that our model 
underestimates the concentration of local subhalos. 
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